Pion wave function from lattice QCD vs. chiral quark models 
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Abstract 

We analyze the equal-time Bethe-Salpeter quark wave function of the pion obtained from a quenched lattice QCD 
calculation with delocalized quark interpolators. We find that the result agrees remarkably well with the predictions 
of the Nambu-Jona-Lasinio model in all channels. We choose the quenched lattice QCD, since it is closer to the large- 
iV c limit of the Nambu-Jona-Lasinio model. We also show how transversity information, relevant for the light-cone 
physics, can be obtained from our equal-time rest-frame lattice calculations. 
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1. Introduction 

In this paper we examine the pion quark wave 
functions in the Nambu-Jona-Lasinio (NJL) model 
and in quenched lattice QCD, and confront both 
results. 

Hadronic wave functions encode important in- 
formation on bound states in strong interaction 
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physics; in particular, they provide the amplitude 
for a composite hadron to have quarks in a given 
momentum state or, equivalently, at a certain space- 
time distance. For systems with heavy quarks non- 
relativistic quantum mechanics applies and particle 
number is conserved, thus much of our understand- 
ing is directly based on wave functions. However, 
as a matter of principle the wave functions cannot 
be directly measured experimentally and one must 
instead resort to form factors, decay widths, or 
momentum distributions. For light quark systems 
particle creation may occur, demanding a field- 
theoretic framework where further complications 
arise; even in the simplest meson case, relativistic 
invariance requires that one uses the conventional 
Bethe-Salpeter amplitude [T] with fixed number of 
quark field operators, a reminiscent of the approx- 
imated parton picture point of view, emphasized 
by the light-cone approaches [21314] . Color gauge 
invariance requires that one additionally includes 
link operators [5 6J. For the pion, the spontaneously 
broken chiral symmetry is a basic dynamical ingrc- 
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dient in the determination of its nonperturbative 
quark structure. It appears via the pertinent axial 
Ward-Takahashi identities |H7j (for a review see, 
e.g., [8] and references therein). These important 
constraints are implemented in relativistic field- 
theoretic chiral quark models, such as the NJL 
model [S]. The regularization needs to be carefully 
handled (for a relevant review see, e.g., [10]). 

On the other hand, lattice QCD solves the bound 
state problem in a fundamental way. It is thus possi- 
ble to make a first-principle non-perturbative deter- 
mination but at the expense of breaking continuum 
symmetries, such as the Lorentz invariance and, 
quite often, chiral symmetry, due to the finite lattice 
spacing. After a pioneering study |llj the hadronic 
wave functions have been analyzed on the lattice on 
a number of occasions [12ll3ll4ll5il6ll7ll8ll9l20] . 
The axial Ward-Takahashi identities can be exactly 
implemented on the discrete Euclidean lattice as 
shown by Ginsparg and Wilson [21] (see Ref. [22|23j 
for a recent practical implementation; here we use 
the same method) , enabling realistically small pion 
masses. 

2. Bethe-Salpeter amplitude 

Since we attempt a comparison between a lattice 
calculation, where only gauge invariant objects are 
defined, and a quark model calculation with no ex- 
plicit mention of the color gauge symmetry, some 
remarks delineating the scope and meaning of such 
a comparison are in order before presenting the ac- 
tual calculations. 

The Bethe-Salpeter vertex (or the quark- 
antiquark wave function) of the pion is given by 

X b q (p) = -i J d 4 xe-^(0\T{q(x)q(0)}\Mq)), (1) 

where q(x) is the spinor field operator carrying fla- 
vor and color, 1 7Tj, (q) ) is the pion state with isospin 
b and on-shcll four- momentum q, q 2 = m~, and p 
denotes the momentum of the quark field after the 
Fourier transform. While chiral quark model calcu- 
lations are naturally formulated in the momentum 
space, the basic objects in the Euclidean lattice cal- 
culations are correlation functions in the coordinate 
space defined in Sect. [4] which are gauge and renor- 
malization group invariant at all Euclidean times. 

The most general form of the quark- antiquark cor- 
relator allowed by the symmetries [24125] has the 
structure 



(Q\T{ q (x)q(0)}Mq)) = ^n (2) 

where the wave functions a = P,A,T, depend 
on the Lorentz-invariant variables x 2 , x ■ q, and q = 
m 2 . The LHS is the inverse Fourier transform of 
the vertex function in the Bethe-Salpeter amplitude 
([T]), which is finite and undergoes cc-independent 
multiplicative renormalization. Thus, the ratios 
\E , (x)/\l/(0) become cut-off independent as the cut- 
off is removed, which on the lattice means the lattice 
spacing a — > 0. In other words, we are studying the 
rcnormalization-group- invariant object. 

The definition (TT]) is satisfactory for chiral quark 
models. In QCD, however, it is appropriate only in 
the fixed point Fock-Schwinger gauge, x^A^x) = 
(the special case being the light-cone gauge 
n^A'Hx) — 0), where the standard derivatives, d^, 
and the covariant derivatives, = + igA^, 
coincide. On the lattice, fixing the gauge has the 
problem of the Gribov copies, as there exists no 
complete gauge fixing for nonabelian nonperturba- 
tive theories. Moreover, Elizur's theorem prevents 
non-vanishing vacuum expectation values of gauge 
variant operators. Non-gauge- invariant bilinear op- 
erators are made gauge invariant by joining them 
with a link operator, however path dependence sets 
inPI Furthermore, gluons carry momentum in the 
pion and any gauge invariant but path-dependent 
definition will yield different results (see Ref. [57] 
for a discussion of various possibilities). For defi- 
niteness, we choose the straight line prescription for 
the path. We also undertake a smearing procedure 
of the link, as described in Sect. [4] 

Another important issue concerns the compari- 
son with either the quenched or dynamical results, 
with the full inclusion of the fermion determinant. 
The NJL Lagrangian models to the one-quark-loop 
level correspond to a large- N c approximation. On 
the other hand, at large N c the fermion determinant 
is suppressed in QCD, explaining why mesons are 
stable in that limit [28129] , The quenched approx- 
imation contains all the leading- iV c and a piece of 
the subleading in N c contributions, which is actually 
suppressed for heavy quarks. That means that pion 
loops are 1 /^-suppressed, although not all of the 



See, e.g., Ref. 26 for an illustration within nonlocal chi- 
ral quark models. The issue reflects the standard operator 
ordering ambiguity between and p M . Only if there is a 
gauge fixing where the link operator becomes the identity, 
the path-independence is guaranteed. 
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1 /N c contributions come from pion loops [3D] . More- 
over, besides including the fermion determinant, one 
should also consider higher Fock states (qq) 2 , etc., 
in the wave function when comparing to the dynam- 
ical lattice results. 

In the present paper we consider the quenched lat- 
tice results and restrict to the NJL model, as the sim- 
plest prototype of a chiral quark model. Our study 
can be viewed as a useful quantitative test of pion 
wave functions which can be carried out for other 
chiral quark models at any level of sophistication. 
We recall here an early comparison of the instanton- 
model hadronic wave functions to lattice results [31] . 

Taking the charged pion for definiteness of the no- 
tation, multiplying ([3]) with appropriate Dirac and 
isospin matrices, and taking the traces yields the re- 
lations directly used in our evaluation: 

(0|J(0)tT6u(s)|7r + (g)) = V2*p, (3) 
(0\d(p)i^u{x)\K + {q)) = vV*A, 
(0\d(0) l5 a^u(x)\TT+(q)) = V2(q"x u - q u x^ T . 

Our choice of the kinematics is q = (m w ,0) and 
x = (0, r). In order to extract the instant-form wave 
functions f „ from the lattice we will consider in 
Sect. 2] the matrix elements 

(0|d(0)r o «(Q, r)|7r+(m OT , 0)> - y/29 a (r), (4) 

with the explicit form of the vertices Tp =175, T a = 
il57o/m n and T T = j 5 a Qi r i /(m n r 2 ). 

3. Wave functions from the NJL model 

A convenient way to determine the pion wave 
function in chiral quark models is by exploiting the 
axial Ward-Takahashi identity (see, e.g., [TU] for the 
details), relating the quark propagator S(j>) and 
the irreducible vertex function, T^ a (p + q,p), cor- 
responding to the axial current. The spontaneous 
breaking of the chiral symmetry generates a con- 
stituent quark mass M given by the gap equation, 
yielding S(p) — i/(j/> — M — to). The pion wave 
function is extracted from the pion pole of T^ a in 
the form of an unamputated vertex function [10] , 



(a) 



x b q ( P ) = 



j) — M — to 



M — m 



(5) 



The pion-quark coupling constant, g^qq, satisfies the 
Goldberger-Treiman relation, g^ qq = M/f^, with f„ 
denoting the pion decay constant. 



q=0, q =m n 



x=r, t=0 
x=0, t=0 



(b) 

q=0, q =m n 



p fixed 



Fig. 1. Fcynman diagrams for the one-quark-loop evaluation 
of the pion wave function at rest (a) in the coordinate space 
and (b) in the momentum space. 

The approach is non-renormalizable and re- 
quires introducing a finite cut-off which should 
fulfill a number of requirements concerning the 
gauge and Lorentz invariance as well as causal- 
ity. This is actually crucial for many applications, 
including the determination of Chiral Effective 
Lagrangeans [32 33 , Parton Distribution Func- 
tions |34|35] , dispersion relations for two point 
functions [36 37], Parton Distribution Amplitudes 
and the light-cone wave functions |38I39I40|4"T] . 
Generalized Parton Distributions 42 43j, Gener- 
alized Form Factors [44], the photon Distribution 
Amplitude 45J, or the pion-photon Transition Dis- 
tribution Amplitude [46|47|48|49|50] . For practical 
calculations the so-called bosonized form is more 
convenient and we refer to 10J for further details. 

We apply the simplest twice-subtracted version of 
the Pauli-Villars regularization. For an observable 
A it amounts to the replacement M 2 — > M 2 + A 2 , 
followed by the subtraction 

A| reg = A(A 2 = 0)-A(A 2 )+A 2 ^l. (6) 

The model has three parameters, which can be 
traded for f„, m^, and M. At a fixed value of M we 
determine m = mo and A by fixing and / w [10] 
to their physical values mjj llys = 139 MeV and f n = 
93 MeV. In this work we use M = 300 MeV, A = 
790 MeV, and too = 8.2 MeV. 

The lattice simulations are performed at m„ > 
m P hys , hence we need to increase accordingly the 
pion mass in the model. At not-too-large values of 
TO,r this may be conveniently achieved via the Gell- 
Mann-Oakes-Renner relation, m\ oc m, from where 
to = TOo(TO 7r /TO,P hys ) 2 . This value of to is actually 
taken for the lattice values of m n . 

Now we evaluate the quantities of our interest, 
^ a (r) = ^ a (0, — r 2 ). The calculation proceeds ac- 
cording to the diagrams of Fig. [T] Standard Feyn- 
man rules yield 



*«(r) 



d 3 p 
(2tt) 3 



dp 



Tr [r a 5 # p^7rgg'75 Sp— q\ 
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(7) 



Then we perform the integration over po and carry 
out the Fourier-Bessel transform over p, with p = 
\p\. The result is 



xT , . _ [dpp 2 . ( 2N c g^p 2 + {M+m) 2 
P[T> ' 2^ 2 MPr > p 2 + (M + m) 2 - m 2 j4 



f dpp 



2tt 2 



g-irqq(M + 771 ) 



^p 2 + {M + m) 2 (p 2 + (M + m) 2 - m 2 / '4) rcg ' 

(8) 

cQirqq 



f dp jp 1 p 



y/p 2 + (M + m) 2 (p 2 + (M + m) 2 - m 2 /4) 

The resulting ratio ^ a (r)/^ a (0), obtained numeri- 
cally from Eqs. ([8|), is plotted in Fig. [2] In the chiral 
limit (m^ = 0, m = 0) one can carry out the inte- 
gration in Eq. ([8]) analytically, which yields 

* F (r)=2# T (r) = %^ (-rA 2 K ( ^ A 2 + M 2 r) 
-2y/A 2 + M 2 Xi(v / A 2 + M 2 r) + 2MK x {Mr 
g^ qq MN c ( A 2 rK 1 (VA 2 + M 2 r) 



4tt 2 



VA 2 + M 2 

-2K ( VA 2 + M 2 r^j + 2K (Mr)) , (9) 

where Kq and K\ are the modified Bessel functions. 
This leads to the following asymptotic behavior at 
r — ► oo: 



-Mr 



r-3/2 



■j "Mr) 



-Mr 



„l/2 



(10) 



We note an exponential fall-off and a longer tail in 
the A channel than in the P and T channels. 

One may also compute the two-dimensional 
Fourier-Bessel transform of Eq. passing from r 
to the transverse momentum kx, which then yields 
the transverse-momentum light-cone wave func- 
tions integrated over a. We find in the chiral limit 



*pO t ) = 2c j7rqq N c 

(. (k 2 +A 2 + M 2 
X (M k 2 + M 2 



(11) 



A 2 



k 2 + A 2 + M 2 



^ A (k 7 



2* T (fc T ) 
M 



2g nqq A i MN c 



{k 2 -, + M 2 )(k\ + A 2 + M 2 ) 



Note that while in the coordinate representation 
\& p(r) ~ $r(f), hi the fcp representation '^(fcp) ~ 

4. Wave function from the lattice 

We now turn to the quenched lattice calculation 
of Hf a {r), defined as (g]) 



y a (r) = l/V2(0\O:\Tr+(q = 0)) } 
O r a (y,t)=d(x,t)V{G}r a u(y + r,t), 



(12) 



and illustrated in Fig. [Ha). "The interpolator anni- 
hilates quark and anti-quark at given time t and 
at distance r apart. The Lorentz invariance ensures 
that ^a^r) is independent of the arguments y and 
t. The gauge link V[G\ from y to y + r ensures the 
gauge invariance and it depends on the choice of the 
path, as discussed in Sect. [H We choose a straight 
path between y and y + r. Since a straight path is 
unique only along the lattice directions r = Naei, 
we evaluate vE f Q (r) only for integer multiples of lat- 
tice spacing a, i.e., r = Na. Then V[G] is a product 
of N gauge links. For a = P, A we use in fact the 
interpolator (1121) averaged over six lattice points at 
the distance r = Na. 

The correlation function is a basic object on the 
Euclidean lattice, which allows one to extract the 
pion mass , as well as its coupling to a given inter- 
polating composite operator (0|O^|tt). For our pur- 
pose we compute a correlation function with a pseu- 
doscalar point source O P =0 at time and a delocal- 



ized sink with a 
time t = —ir: 

C r a (T) = [ d 3 ye 



P, A, or T, O r a at the imaginary 



a l iv / 



0:(r,y)O r P ^ (0,0)|0) (13) 



where we project on the total momentum q = 0. 

The quantities related to the pion can be ex- 
tracted after the complete set of physical states is 
inserted into (113|) . yielding 



\O r M\n q )(n q \O r p = a \m r - En(Q) r 



2E n (q) 



(14) 



The ground-state pion at rest dominates at large r, 
C-(r) r -^?< e -"^, (15) 
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\OZ\ir) 



WOp =ot |o) 



2to w 

From a single exponential fit to the lattice correla- 
tion functions at large r we extract m n and w r a for a 
range r = [0, a, .., 4a]. We use several current quark 
masses to, which corresponds to m v in the range 
345 - 740 MeV. 

The extracted pion mass does not numerically de- 
pend on r nor the choice of the channel P, A, or T, 
as expected. The ratios ^> a (r)/ fy a (0) for P and A 
channels are extracted from the identity 



*„(r) 



\O r = 



(16) 



*a(0) 

and plotted in Fig. [2 The wave function can not 
be evaluated at r = on the lattice and we normalize 
it arbitrarily, such that the NJL model results and 
lattice values at r = a agree (c/. Fig. [2]). 

Finally, we provide some details of the lattice sim- 
ulation. We use quenched lattice QCD, since it is 
closer to the large- N c limit of the NJL model, as ex- 
plained in Sect.dl We use 100 gauge configurations 
generated by the Liischer-Weisz gauge action [ST], 
as described in [52] . 

We use so called HYP smearing on the gauge con- 
figurations [53] . This replaces a gauge link between 
the neighboring points on the lattice with a "fat" 
link. The fat link is a sum of links within hyper- 
cubes attached to the original link only. Since this 
smearing is local, it does not destroy short distance 
quantities, while it smooths out the large local fluc- 
tuations of the gauge field. 

The lattice spacing a ~ 0.148 fm is determined 
from the Sommer parameter, and the lattice volume 
is 16 3 x 32. We employ the chirally improved valence 
quarks 22 23 52J, which have good chiral properties. 

5. NJL vs. quenched lattice 

Our results are presented in Fig. [2] with the lines 
showing the model and the points the lattice calcu- 
lation. In the P and A channels we normalize $F a (r ) 
to its value at the origin, while for the T channel we 
normalize to the model prediction at r = a. We note 
that the results of the NJL model agree remarkably 
well with the lattice determination in all three chan- 
nels. The dependence on m ff is very weak. The insen- 
sitivity of the lattice wave functions to the value of 
the pion mass confirms previous findings [14ll5llt?] . 
As expected from Eq. ([9]) , the pseudoscalar and ten- 
sor wave functions are identical within the error 
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Fig. 2. The components of the rest-frame pion wave function, 
normalized to unity at the origin, ^/ a (r)/^!a(0), for S, A, and 
T channels, evaluated in the NJL model at M = 300 MeV 
and compared to the quenched lattice data. The shown points 
are at = 345 and 575 MeV, while the model calculation 
includes also the case of the physical pion mass. 

bars, while the axial wave function is found to be 
wider. This agrees with the different asymptotic be- 
havior of Eq. (TIT)]) . 

We have also checked that the model results are 
almost independent of M when this parameter is in 
the reasonable range, M = 250 — 350 MeV. 



6. Rest-frame kinematics and transversity 

We end this paper with general remarks relat- 
ing the equal-time wave functions presented above 
to the light-cone wave functions in the impact- 
parameter representation, also known as transver- 



5 



sity wave functions. Formally, we can write the 
transformation 

1 

(Q\T {q(x)q(0)} \ir b (q)) = % J^J dae ^-^ x 

o 

[-<P P {a, x 2 ) + A (a, x 2 ) - ia^q^x u ^ T (a, x 2 )], 

(17) 

where a is the Feynman parameter [24]. The pres- 
ence of a gauge link operator between the quarks is 
implicit. Comparison to ([3]) yields by construction 

l 

V a (x-q,x 2 ) = I dae^ 2a -^ x <i> a (a : x 2 ). (18) 
o 

As a matter of principle, all scalars \& a in Eq. (|17p 
depend on the scalar variables x 2 , q 2 , and x ■ q, 
hence we are free to choose any reference frame to 
evaluate ^ a . In the rest-frame, or equal-time (ET) 
kinematics, used in the previous sections of the pa- 
per, we take x — (0, r) and q = (m T ,0), whence 
x 2 = —r 2 and x ■ q — 0. In the infinitc-momcntum- 
frame kinematics (q , q) = lim g ^ 00 (-\/m^ + q 2 , q z ) 
and on the light cone (LC), where x + = 0, one 
has x ■ q = q + x~, x 2 = — r 2 . The parameter a 
acquires the meaning of the light-cone momentum 
fraction of pion carried by one of the quarks. By com- 
paring the two calculations we find ^ T (0, — r 2 ) = 
/ 1 dae t ( 2Q - 1 )« +x_ $L c (a,-r 2 ).Forthe chosen kine- 
matics q + x~ = q • x = 0, hence 

l 

^ T (0,-r 2 ) = J da^ a c ( ai -r 2 ). (19) 
o 

In other words, our rest frame calculation allows for 
a direct determination of the transverse-coordinate 
(impact parameter) dependence of the light-cone 
wave-function integrated over the a parameter. A 
similar property for the Generalized Parton Dis- 
tributions was suggested for the nucleon [55] and 
the pion [53]. The relation between the equal-time 
and light-front wave functions has been analyzed 
recently |55j in the momentum space, where the 
transversity relation (|19[) cannot be explicitly seen. 

It would also be useful to verify the equal time- 
light cone transversity connection on the lattice. 
While there exist transverse lattice calculations [55] , 
their focus was on the Distribution Amplitude, 
ty(a,0) = </?(a), leaving out the transverse depen- 
dence. Some results where also presented in [57] . 



The transversity relation (ITS]) is explicitly verified 
for the N JL light-cone wave function [39 ■ 

7. Conclusions 

Here are our main points: 

- The leading- N c chiral quark model interpretation 
of the quenched lattice data is not only qualita- 
tively correct, but also remarkably accurate. This 
is yet another manifestation of the fact that the 
spontaneously broken chiral symmetry is the key 
dynamical factor in the pion dynamics. The qual- 
ity of the agreement suggests that the 1 /N c con- 
tributions in the quenched smeared lattice simu- 
lations are small. 

- The dependence of both the lattice and model 
results on the value of the pion mass is very small 
in the broad tested range m w = 345 — 740 MeV. 

- The pseudoscalar and tensor quenched pion wave 
functions are equal on the lattice within the error 
bars. In the model they are equal in the chiral 
limit and nearly equal for the used values of m n . 

- The asymptotic fall-off of the pion wave functions 
in the model is exponential, ~ exp(— Mr)/r p , 
with p = 3/2 in the pseudoscalar and tensor 
channels, while p = 1/2 in the axial channel. This 
qualitatively complies to the lattice data, where 
the axial wave function exhibits a longer tail. 

- A general result, which originates from the 
Lorentz invariance, concerns the utility of the 
equal-time rest-frame smeared lattice simulations 
to determine the transversity information relevant 
for the light-cone physics. The integrated light- 
cone (infinite-momentum) pion wave functions 
in the impact-parameter space, J Q da$>^ c (a, 6), 
coincide with our equal-time rest-frame wave 
functions evaluated at the same quark-antiquark 
separation, ^ T (0, —r 2 ) \ r=b . 

One of us (SP) thanks Bojan Golli, Christian B. 
Lang, and Mitja Rosina for useful discussions. We 
are grateful to the Bern-Graz-Regensburg Collab- 
oration for providing the quenched lattice config- 
urations and the quark propagators, indispensable 
for this work. They were generated in the Leibniz 
Rechenzentrum, Garching and the ZID, Graz. 
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